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Abstract 

We develop a "variational mass" expansion approach, recently introduced in the 
Gross-Neveu model, to evaluate some of the order parameters of chiral symmetry 
breakdown in QCD. The method relies on a reorganization of the usual perturbation 
theory with the addition of an arbitrary quark mass m, whose non-perturbative be- 
haviour is inferred partly from renormalization group properties, and from analytic 
continuation in m properties. The resulting ansatz can be optimized, and in the chi- 
ral limit m — > we estimate the dynamical contribution to the "constituent" masses 
of the light quarks M^^s] the pion decay constant and the quark condensate 
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1 Introduction 



A still challenging question in strong interaction physics is the derivation of the low- 
energy properties of the spectrum from "QCD first principle" , due to our limited present 
skill with non-perturbative physics. At very low energy, where the ordinary perturbation 
theory cannot be applied, Chiral Perturbation Theory and Extended Nambu-Jona- 
Lasinio models 0, Q give a consistent framework in terms of a set of parameters that 
have to be fixed from the data; yet the bridge between those effective parameters and 
the basic QCD degrees of freedom remains largely unsolved. Although lattice QCD sim- 
ulations recently made definite progress in that direction, the consistent treatment of 
dynamical unquenched quarks and the chiral symmetry remains a serious problem. 



In this paper, we investigate a new, semi-analytical method, to explore how far the 
basic QCD Lagrangian can provide, in a self-consistent way, non-zero dynamical quark 
masses, quark condensates, and pion decay constant, in the limit of vanishing Lagrangian 
(current) quark masses. Such a qualitative picture of chiral symmetry breakdown (CSB) 
can be made more quantitative by applying a new "variational mass" approach, recently 
developed within the framework of the anharmonic oscillator , and in the Gross-Neveu 
(GN) model [0, ||. The starting point is very similar to the ideas developed a long time 
ago and implemented in various different forms in refs.0, [|10]. There, it was advocated 
that the convergence of conventional perturbation theory may be improved by a varia- 
tional procedure in which the separation of the action into "free" and "interaction" parts 
is made to depend on some set of auxiliary parameters. The results obtained by expand- 
ing to finite order in this redefined perturbation series are optimal in regions of the space 
of auxiliary parameters where they are least sensitive to these parameters. Recently there 
appeared strong evidence that this optimized perturbation theory may indeed lead to a 



rigorously convergent series of approximations even in strong coupling cases |TT 



An essential novelty, however, in and the present paper, is that our construc- 

tion combines in a specific manner the renormalization group (RG) invariance with the 
properties of an analytically continued, arbitrary mass parameter m. This, at least in a 
certain approximation to be motivated, allows us to reach infinite order of the variational- 
perturbative expansion, and therefore presumably optimal, provided it converges. Our 
main results are a set of non-perturbative ansatzs for the relevant CSB quantities, as func- 
tions of the variational mass m, which can be studied for extrema and optimized. Quite 
essentially, our construction also provides a simple and consistent treatment of the renor- 
malization, reconciling the variational approach with the inherent infinities of quantum 
field theory and the RG properties. 

Before proceeding, let us note that there exists a quite radically different attitude 
towards CSB in QCD, advocating that the responsible mechanism is most probably the 
non-perturbative effects due to the instanton vacuum |T2|, or even more directly related 
to confinement [|13|. However, even if the instanton picture of CSB is on general grounds 
well motivated, and many fruitful ideas have been developed in that context^], as far 
as we are aware there is at present no sufficiently rigorous or compelling evidence for 

^ For a recent complete review see Q . 

^See e.g ref. [Ol for a review and original references. 
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it, derived from "first principle". In any event, it is certainly of interest to investigate 
quantitatively the "non-instantonic" contribution to CSB, and we hope that our method 
is a more consistent step in that direction. 

2 Dynamical quark masses 

In what follows we only consider the SU{nf)L x SU{nf)R part of the chiral symmetry, 
realized by the QCD Lagrangian in the absence of quark mass terms, and for nj = 2 or 
n/ = 3 as physically relevant applications. 

Following the treatment of the anharmonic oscillator and its generalization to the 
GN model 0, |[, let us consider the following modification of the usual QCD Lagrangian, 

Uf Uf 

Lqcd L^QCDi.9Q = 0, mo = 0) - rriQ^^q^qi + Lgcnigl ^ xg^) + x mo^^Mi , (1) 

i=l i=l 

where L^qcd designates the ordinary QCD interaction terms, and a; is a convenient "in- 
terpolating" expansion parameter. This formally is equivalent to substituting everywhere 
in the bare Lagrangian, 

mo^mo(l-x); gl-^glx, (2) 

and therefore in any perturbative (bare) quantity as well, calculated in terms of mo 
and g^. Since the original massless QCD Lagrangian is recovered for x — 1, mo is to 
be considered as an arbitrary mass parameter after substitution (Q). One expects to 
optimize physical quantities with respect to mo at different, possibly arbitrary orders of 
the expansion parameter x, eventually approaching a stable limit, i.e flattest with respect 
to mo, at sufficiently high order in x. 

However, before accessing any physical quantity of interest for such an optimization, 
the theory should be renormalized, and there is an unavoidable mismatch between the 
expansion in x, as introduced above, and the ordinary perturbative expansion as dictated 
by the mass and coupling counterterms. Moreover, it is easy to see that at any finite order 
in the x expansion, one always recovers a trivial result in the limit m — > (equivalently 
a; — *^ 1), which is the limit in which to identify non-zero order parameters of CSB. 

These problems can be circumvented by advocating a specific ansatz, which resums 
the (reorganized) perturbation series in x and is such that the limit x — 1 no longer gives 
a trivial zero mass gap. As was shown in detail in ref. [§], the ansatz for the dynamical 
mass is most easily derived by following the stepsQ: 

i) Consider first the general solution for the running mass, given as 

m(/.')=m(/z) «-p{-/J7^^w} 

in terms of the effective coupling g{ii), whose RG evolution is given as fidg{fi)/dfi = l3{g), 
and 'Jmig) = — (/i/m)(i(m(/i))/(i/i. Solving (||) imposing the "fixed point" boundary 
condition: 

M = m{M), (4) 
^See also ref. for a detailed derivation in the QCD context. 
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at two- loop RG order we obtain, after some algebra (we use the normalization f3{g) 

-hoQ^ - hiQ^ , ^rn{9) 



M2=m /"W r 



1 + ^^' 



JO TO 

26i 2i)o 



(5) 



where m = rn^jj), g = g{ii) {ji = ji^f^e "t^/^), and / = g^ j g^{M-i) satisfies 



.2,„^"2 , &i.2,„riii^ 



/ = l + 26o^Mn^ + -i^Mn 



/ 



(6) 



(note in (|^) and (P) the recursivity in both / and M^). The necessary non-logarithmic 
perturbative corrections to those pure RG results are then consistently included as 

A/f.M.(l + |7„^+l^|; + 0(/)), (7) 



where the complicated two-loop coefficient K was calculated exactly in ref. |TB[. Equa- 
tion. (^ defines the (infrared-convergent, gauge-invariant) pole mass [|T7| M^", in terms 

to resum the leading (LL) and 



of the MS mass at two-loop order, and can be shown 
next-to-leading logarithmic (NLL) dependence in m to all orders. 

ii) Perform in expressions (^, (|]), (0) the substitution fn fnv, and integrate the re 
suiting expression, denoted by M|'(t>), according to 

1 f dv 



im J V 

where the contour is around the negative real v axis. 



(8) 



In [H it was shown that the previous steps correspond (up to a specific renormalization 
scheme (RS) change, allowed on general grounds from RG properties) to a resummation of 
the X series as generated from the substitution (0)0. Moreover this is in fact the only way 
of rendering compatible the above x expansion and the ordinary perturbative one, thus 
obtaining finite results. Actually the resummation coincides with the exact result in the 
large- limit of the GN model. Now, since the summation can be formally extended to 
arbitrary RG orders |^, including consistently as many arbitrary perturbative correction 
terms as known in a given theory, in the QCD case we make the assumption that it 
gives an adequate "trial ansatz", to be subsequently optimized in a way to be specified 
next. After appropriate rescaling of the basic parameters, g and m, by introducing the 
RG-invariant basic scale A^g^ (at two- loop order), and the convenient scale-invariant 
dimensionless "mass" parameter 



2^ [2h^-g 



-2l-7o/(2feo) 



1 + '^f 

bo 



7o/(2feo)-7i/{2fei) 



we end up with the following dynamical mass ansatz: 



m ) 



-c 



A- 



MS 



2m 



dy 



y/m 



F^iy)[C + Fiy)] 



B 



1 + 



Ml 



+ 



F{y) F^{y)^ 



(9) 



(10) 



V is related to the original expansion parameter a: as a; = 1 — v/q, q being the order of the expansion. 
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where y = m"f , and F is defined as 

F{y) ^ \n[y] - A ln[F(|/)] - {B - C) \n[C + F{y)l (11) 

with A = 7i/(26i), B = 70/(260) - 7i/(26i), C = bi/{2bl), in terms of the RG coefficients 
|T^. Finally the perturbative corrections in (p!0| ) are simply given as A^i = (2/3) (70/260) 
and = K/{2boy. 

Observe in fact that, were we in a simplified QCD world, where there would be no 
non-logarithmic perturbative contributions (i.e. such that A4i = M.2 = ■ ■ ■ = in (p!0|)), 
the latter ansatz would then resums exactly the x variational expansion. In that case, ([To|) 
would have a very simple behaviour near the origin m" — > 0. Indeed, it is easy to see that 
(|ll|) admits an expansion F{y) ^ (7(^^'^)/"^ y^/"^ for ?/ ^ 0, which immediately implies 
that (|l^) would give a simple pole at ?/ ^ 0, with a residue giving M2 = {2C)~^ ^Ws- 
Moreover one can always choose an appropriate renormalization scheme in which 62 and 72 
are set to zero, as well as all higher order coefficients, so that there are no other corrections 
to the simple above relation. 

Now, in the realistic world, A^i, M.21 etc can presumably not be neglected. We can 
nevertheless expand (0) near m" for any known non-zero A^i, using 

— I dyey/^' y'^ = ^"^r ^'T , (12) 
2m J ^ ^ T[-a\ ' ^ ' 

and the resulting Laurent expansion in (m")^/"^ may be analysed for extrema and opti- 
mized at different, in principle arbitrary {m"Y/^ orders. An important point, however, 
is that the perturbative corrections do depend on the RS choice, as is well known. Since 
the pure RG behaviour in ([TOD already gives the order of magnitude, M ~ const x A;|g^, 
we can hope that a perturbative but optimized treatment of the remaining corrections is 
justifed. In other words we shall perform an "optimized perturbation" with respect to 
m" around the non-trivial fixed point of the RG solution. 



To take into account this RS freedom, we first introduce in (|TD|) an arbitrary scale 
parameter a, from ft ^ a ft. Accordingly the perturbative coefficients Aii in (lU) take a 
logarithmic dependence in a, simply fixed order by order from (H)-(|3) and the require- 
ment that ([ToD differs from the original MS expression only by higher order terms. The 
a-dependence will eventually exhibit a non-trivial extrema structure and we shall also 
optimize the result with respect to oQ. Actually there are other possible changes of renor- 
malization prescriptions affecting expression (p!OD in addition to the a dependence, which 
may be taken into account as well. More precisely, the second coefficient of ■ym{g), 7i, do 
depend on the RS choice, even in MS schemes [0. As it turns out, this additional RS 
freedom is very welcome in our case: in fact, the previous picture is invalidated, due to 
the occurence of extra branch cuts in the y plane at Re[ycut] > 0, as given by the zeros 
of dy/dF from (pH]) (in addition to the original cut on the negative real y axis). This 
prevents using the expansion near the origin, eq. (|12D, since it would lead to ambiguities 



^This procedure indeed gave very good results [g in the GN model, where in particular for low values 
of N the optimal values found, Oopt, are quite different from 1. 



4 



of (9(exp(i?e[?/]/m")) for m" OQ. The specific contour around the negative real axis was 
suggested by the known properties of the large limit of the GN model, and it is not 
surprising if the analytic structure is more complicated in QCD. However, the nice point 
is that the actual position of those cuts do depend on the RS, via A (71) in (]Tl|). Defining 
7^ = 7i + A71, we can choose Re[ycut] — for A71 ~ 0.00437 (0.00267) for Uf = 2 {uf = 



3), respectively. We therefore consider |T5[ the general RS change 

' m (1 + B^f + B,f) ■ = f {1 + A,f + A,f) 



m 



(13) 



(implying A71 = 26o-Bi — 7o^i); and optimize with respect to this new arbitrariness^]. 

However one soon realizes that our extension of the "principle of minimal sensitivity" 
(PMS) defines a rather complicated optimization problem. Fortunately, we can study 
this problem within some approximations, which we believe are legitimate. Since the 
ansatz ([To|) (with the above RS change understood, to make it consistent) would indeed 
be optimal with respect to m for vanishing perturbative non-logarithmic corrections, 
Aii = 0, we shall assume that the expansion for small m" is as close as possible to an op- 
timum, and define the m" limit by some relatively crude but standard approximation, 
avoiding numerical optimization with respect to m". The approximation we are looking 



for is not unique: given (10), one could construct different approximations leading to a 
finite limit for m" ||^ . Here we shall only demonstrate the feasibility of our program 
in the simplest possible realization. In fact, since we shall anyhow optimize with respect 
to the RS dependence we assume that it largely takes into account this non-uniqueness 
due to higher order uncertainties. 

Fade approximants are known to greatly improve perturbative results [^] and often 
have the effect of smoothing the RS dependence. We thus take a simple Fade approximant 
which by construction restitutes a simple pole for F — (i.e m" 0) in (0), and gives 



M 



Fade , 



a, A7i,5i,m" ^ 0) = (2C)-^ a exp{^} 

/On 



1 - 



A^2(a,A7i,5i) 



(14) 



We have performed a rather systematic study of the possible extrema of ([14|) for arbitrary 
a, Bi (with A71 fixed so that the extra cuts start at Re[y] ~ 0). We obtain the flattest 
such extrema for a ~ 2.1, Bi ^ 0.1, which leads to the result 



M„^r(m" 



0) ^ 2.97 A^(2) 



(15) 



for Uf = 2. Similarly, we obtain M„^f<^(m" ^ 0) ^ 2.85Ajg5(3) for Uf = 3. Note that 
these values of the dynamical quark masses, if they are to be consistent with the expected 
range |[ of Mdyn ^ 300-400 GeV, call for relatively low A^ ~ 100-150 GeV, which is 
indeed supported by our results in the next section. 



^ The origin of those singularities is rather similar to the ambiguities related to renormalons . An 
essential difference, however, is that the present singularities occur in the analytic continuation of a mass 
parameter rather than a coupling constant, and it is possible to move those singularities away by an 
appropriate RS change, as we discuss next. See ref. for an extended discussion. 



We also impose a further RS choice, 62 = Oi 72 = 0; which fixes A2, B2 in (|13| ) and guarantees 
that our two-loop convention for Aj^g remains unaffected. Note, however, that (|13|) implies A^g^ — > 
AjTg exp{^} = A'. In what follows we express the results in terms of the original A 



MS- 
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3 Composite operators and F-j^ 

We shall now generalize the ansatz (0) for the pion decay constant F^^. The main idea 
is to do perturbation theory around the same RG evolution solution with the non-trivial 
fixed point, as specified by the function F in (p^, with perturbative correction terms 
obviously specific to F^. A definition of suiting all our purposes is ^ ^ 

t J d'qe'^-^0\T A;(x)A^(0)|0) = + 0{p,p.) (16) 

where the axial vector current = (^7/^75A*g)/2 (the A*'s are Gell-Mann SU{3) matrices 
or Pauli matrices for n/ = 3, n/ = 2, respectively). Note that according to (|16D, is to 
be considered as an order parameter of CSB ||23| . 



The perturbative expansion of (|T6D for m 7^ is available to the three-loop order, as it 
can be easily extracted from the very similar contributions to the electroweak p-parameter. 



calculated at two loops in |25| and three loops in [26 



The appropriate generalization of (|iy) for F^. now takes the form 



F^ 2"^'"fl^ r dv , " 1 



MS 



<5,(l + ^^ + ^ (17) 

in terms of F{y) defined by eq. (jl^) and where 5,^, «7r(l) and /?7r(l), whose comphcated 




expressions will be given elsewhere are fixed by matching the perturbative MS ex- 
pansion in a way to be specified next. Formula ( [T7| ) necessitates some comments: apart 
from the obvious changes in the powers of F, y, etc, dictated by dimensional analysis, 
note that the perturbative expansion of the (composite operator) {A^Ay) in ([TBI) starts 
at one- loop, but zero order. This leads to the extra 26oF factor in ([17|) , corresponding 
to an expansion starting at (9(l/(?^)0. Another difference is that the perturbative expan- 
sion of (p!6[) is ambiguous due to remaining divergences after mass and coupling constant 
renormalization. Accordingly it necessitates additional subtractions which, within our 
framework, are nothing but the usual renormalization procedure for a composite oper- 
ator, which is (perturbatively) well-defined [|^. The only consequence is that, after a 
consistent treatment of the subtracted terms (i.e respecting RG invariance), the unam- 
biguous determination of the 1/F" perturbative terms in (|T7|) necessitates the knowledge 
of the in + 1) order of the ordinary perturbative expansion. The nice thing, however, is 
that the subtracted terms only affect the values of and but not the jorm of the 



ansatz ([T7|), as soon as the order of the variational-perturbative expansion is larger than 
1 . The consistency of our formalism is checked by noting that the re-expansion of (P!7| ) 
do reproduce correctly the LL and NLL dependence in fn of the perturbative expansion 
of the composite operator to all orders. 

The analyticity range with respect to A71, discussed in section 2, remains valid for 
(|r^) as well, since the branch cuts are determined by the very same relation (|TT]). We can 
thus proceed to a numerical optimization with respect to the RS dependence, along the 

* The 0(1/3^) first-order term cancels anyhow after the necessary subtraction discussed below. 
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same line as the mass case in section 2. Using an appropriate Pade approximant form to 
define the F — > (m" 0) hmit, we obtain the optimal values as 

<ot'(^" - 0) ^ 0.55 A^(2) (0.59 A^(3) ) , (18) 

for nf = 2 (3). With ~ 92 MeV, this gives ^ 157 (168) MeV, for Uf = 3 (2). 



4 (qq) ansatz 



As is well known [pO| , ^, (qq) is not RG-invariant, while m{qq) is; this is thus the rele- 
vant quantity to consider for applying our RG-invariant construction. A straightforward 
generalization of the derivation in section 3 leads to the ansatz 



Agr - (2^°)^/ J {FYA-^[C + Fr^,.) + ) (19) 

where again the coefficients and a(qg)(l) are obtained from matching the ordinary 
perturbative expansion after a subtraction, and will be given explicitely elsewhere |T^ . 



The perturbative expansion, known up to two-loop order p7| , |T5| implies that one only 
knows unambiguously the first order correction 0{1/F) in (|19|), as previously discussed. 
Apart from that, (|T9|) has all the expected properties (RG invariance, resumming LL and 
NLL dependence etc), but a clear inconvenience is that (qq) cannot be directly accessed, 
being screened by tiny explicit symmetry breaking effects due to m 7^ 0. This is of course 
a well-known problem, not specific to our construction. 

However, it is not clear how to consistently include explicit symmetry breaking effects 
within our framework. As amply discussed, in (^) m" is an arbitrary parameter, destined 
to reach its chiral limit m" — >■ 0. Accordingly, m — > for m" — > 0, so that one presumably 
expects only to recover a trivial result, rh{qq) for m" 0. This is actually the 
case: although (|l^) potentially gives a non-trivial result in the chiral limit, namely the 
simple pole residue (= 26o(2C)~'" 6(^gg) a(^gq){a), upon neglecting unknown higher-order 
purely perturbative corrections), when we require extrema of this expression with respect 
to RS changes, using for the m" —>■ limit a Pade approximant similar to the one for 
Ftt, we do not find non-zero extrema. Such a result is not conclusive regarding the ac- 
tual value of {qq){fi), but it may be considered a consistency cross-check of our formalism. 



On the other hand, we should mention that our basic expression (^) does possess non- 
trivial extrema for some m'^p^ 7^ 0. These we however refrain from interpreating in a more 
quantitative way since, within our framework, we cannot give to m" x the meaning 
of a true, explicit quark mass (whose input we in principle need in order to extract a {qq) 
value from (0)). At least, it strongly indicates that it should be possible to extract {qq) 
in the chiral limit, by introducing in a consistent way a small explicit symmetry-breaking 
mass, —mo^exp(li(li, to the basic Lagrangian (P. 



5 Summary 

In this paper we have shown that the variational expansion in arbitrary m", as developed 
in the context of the GN model f^, can be formally extended to the QCD case, apart 
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from the complication due to the presence of extra singularities, which can be however 
removed by appropriate RS change. As a result we obtain in the chiral limit non-trivial 
relationships between A-jg-g and the dynamical masses and order parameters, F^^, rh{qq). 
The resulting expressions in a generalized RS have been numerically optimized, using 
a well-motivated Pade approximant form, due to the complexity of the full optimization 
problem. The optimal values obtained for Mg and are quite encouraging, while for (qq) 
they are quantitatively not conclusive, due to the inherent screening by an explicit mass 
term of this quantity, in the limit m — 0. A possible extension to include consistently 
explicit breaking mass terms in our formalism is explored in ref. [^. 
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